Towards the 4th post-Newtonian Hamiltonian for two-point-mass systems 
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The article presents the conservative dynamics of gravitationally interacting two-point-mass sys- 
tems up to the eight order in the inverse power of the velocity of light, i.e. 4th post-Newtonian 
(4PN) order, and up to quadratic order in Newton's gravitational constant. Additionally, all loga- 
rithmic terms at the 4PN order are given as well as terms describing the test-mass limit. With the 
aid of the Poincare algebra additional terms are obtained. The dynamics is presented in form of 
an autonomous Hamiltonian derived within the formalism of Arnowitt, Deser and Misner. Out of 
the 57 different terms of the 4PN Hamiltonian in the center-of-mass frame, the coefficients of 45 of 
them are derived. Reduction of the obtained results to circular orbits is performed resulting in the 
4PN-accurate formula for energy expressed in terms of angular frequency in which two coefficients 
are obtained for the first time. 
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I. INTRODUCTION 

Calculations of post-Newtonian (PN) expressions 
within general relativity have proved very useful because 
of their full-fledged analytic structures through very high 
orders in the approximation scheme, wherein the nth PN 
order denotes the order (1/c 2 )™ with c denoting the speed 
of light. For compact binaries without spinning compo- 
nents the dynamics is explicitly known through 3.5PN 
order whereby the orders 2.5PN and 3.5PN are of dissipa- 
tive type, see e.g. [H-0|- Partial knowledge is available at 
higher orders. To all orders in 1/c, the n-body dynamics 
is known through linear order in Newton's gravitational 
constant G |8| , and in the order 1/c 8 the terms quadratic 
in G were presented quite recently Q. The terms at the 
4PN level containing logarithms are known too [Iol - [l2^ . 

In case of spinning binaries the conservative dynamics 
is known through 3PN order in both spin-orbit and spinl- 
spin2 couplings or, if the spins are counted 1/c, 3.5PN 
and 4PN order, respectively, see e.g. p^Hl5j . For the dis- 
sipative spin-orbit and spinl-spin2 couplings the highest 
achieved order is 3.5PN (counting spin as of order zero) 

pun. 

Those calculations are very important for future gravi- 
tational wave astronomy because they allow the determi- 
nation of the binary insipral process with high precision, 
[20| . Even last stable binary orbits can be deduced from 
them with high confidence though they are located in 
the limiting regime of PN approximations. On the other 
side, the application of PN results in the effective one- 
body approach (EOB) even allows the transition beyond 
the last stable orbit regime, see e.g. [2l|,[22j]. No doubt, 



the knowledge of any part of the 4PN point-mass dynam- 
ics is of great importance. 

We employ the following notations: x = (a^ 1 ), i = 
1,2,3, denotes a point in the 3-dimensional Euclidean 
space R 3 endowed with a standard Euclidean metric and 
a scalar product (denoted by a dot). Letters a, b are body 
labels (a, b = 1, 2), so x a € M 3 denotes the position of the 
ath point mass. We also define r i2 = x 2 — xi, r i2 = |r i2 1 , 
n = ri 2 /ri 2 ; | • | stands here for the Euclidean length of a 
vector. The linear momentum vector of the ath body is 
denoted by p a = (p a i), and m a denotes its mass param- 
eter. An overdot means the total time derivative. We 
abbreviate 5(x — x a ) by S a . Throughout the paper we 
extensively used the computer-algebra system Mathe- 

MATICA. 



II. THE ADM CANONICAL APPROACH 

In this Section we use units in which 1%-kG = c = 1. In 
the Arnowitt, Deser, and Misner (ADM) canonical ap- 
proach to general relativity [23[ , the constraint equations 
are the crucial equations to be solved for the description 
of the dynamics. Written for many-point-mass systems 
they read 

9 ; 2 (gR • \{m^t "^") 

= (^PaiPaj + m T) 1/2 S a, 
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(2.1b) 
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Here, the 3-metric reads and its inverse is denoted 
by 7 1,7 , g is the determinant of the 3-metric, R is the 
curvature scalar of the time-equal-constant slices, and 



2 



7r lJ the canonical conjugate to gij. The 3-dimensional 
covariant derivative is denoted by | . We have employed 
the ADMTT coordinate conditions, 



9ij 
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1 



6ij + h l3 



TT 



0. 



(2.2) 



where all time derivatives of x a and p a are eliminated 
through lower-order Hamilton equations of motion [24j . 
what is equivalent to performing a canonical transforma- 
tion. For computation of the 4PN matter Hamiltonian 
one needs to use Newtonian and 1PN equations of mo- 
tion. 



where hj,^ is a transverse-traceless (TT) quantity. The 
field momentum 7r y is splitted into its longitudinal Tr^ 
and TT 7r^? T parts, 7r y = -Kij + tt^ t . By solving the 
constraint equations (|2.ip . the ADM Hamiltonian can be 
put into its reduced form, 

ffrcd [x Q , p a , flj?, 7T^ T ] = - J S X A(j). (2.3) 

This Hamiltonian describes the evolution of the matter 
(x a ,p a ) and independent gravitational field (/i^ T ,7r^ T ) 

variables. For the latter, 7r^? T is the canonical conjugate 

An autonomous conservative Hamiltonian can be ob- 
tained through the transition to a Routhian description, 



R 



x a i Pa i ">ij , hij 



= H led - /A4 T ^ T . (2.4) 



TT ij 

Then the matter Hamiltonian reads 

H(x a ,Pa) = R X Q , p a , /l^ T (x a , p a ), ^ T (x a , Pa) 



(2.5) 



III. RESULTS 

We have iteratively solved, by a PN expansion up to 
the 4PN order, the constraint equations (|2.ip for the 
functions <fi and 7^- (the iterative solution of the con- 
straints up to the 3PN order is given in detail in Ref. 
[25|). We have thus obtained the 4PN- accurate Hamil- 
tonian density. The 3-dimensional integral over it we 
have regularized by means of procedures described in 
Appendix B of [25J. To diminish number of terms we 
display here the Hamiltonian at the 4PN order only 
and in the center-of-mass frame defined by the condi- 
tion pi + P2 = 0. We introduce the following reduced 
variables: r = Y\2/(GM) (with r = |r| and n = r/r), 
p = Pi/V, where M = mi + rri2 is the total mass of the 
system and u = mi 7712/ 'M is its reduced mass. We also 
introduce the reduced Hamiltonian H = (H — Mc 2 )//i 
which depends on masses only through the symmetric 
mas ratio v = fi/M (0 < v < 1/4; ^ = is the test- 
mass limit and v = 1/4 holds for equal masses). The 
4PN reduced Hamiltonian reads [here (p 2 ) = p • p and 
(np) = n p] 
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wherein the non computed terms C 4 i(r,p) and C 2 i(r,p) have the structure 

C 4l (r,p) = c 4l i{p 2 ) 2 + c il2 {np) 2 p 2 + c il3 (np) 4 , 
C 2l (r,p) = c 2ll {p 2 ) + c 2i2 (np) 2 , i = 1,2, 



= 1,2, 



(3.1) 

(3.2a) 
(3.2b) 



where c 4 a, c 4i2 , c 4i2 , c 4a3 , c 2 ;2 (i = 1,2) together with 
coi and C02 are not yet computed numerical coefficients. 
The Hamiltonian (|3.1I) is composed of 57 terms, coef- 
ficients of 45 of them are computed. All terms of the 
order 10, 8, and 6 in momenta are calculated. Regu- 
larization of integrals related with these terms does not 
lead to any ambiguity. We have checked that the terms 
of the order 8 in momenta (they are linear in Newton's 
gravitational constant G) are, up to adding a total time 
derivative, compatible with the 4PN Hamiltonian which 
follow from the exact post-Minkowskian Hamiltonian de- 
rived in Ref. Q. We have computed all terms for v = 
and have checked that they agree with the 4PN approxi- 
mation of the Hamiltonian of a test body orbiting around 
a Schwarzschild black hole [see, e.g., Eq. (77) in [25j|]. 
Linear and quadratic in Newton's gravitational constant 
G terms in Lagrangian form were recently derived in Q . 

We have computed all the logarithmic terms of the 
Hamiltonian resulting from divergences of the instanta- 
neous near-zone metric when going to large distances. In 
these logarithmic terms s is a non- fixed constant regular- 
ization scale [more precisley s is the reduced regulariza- 
tion scale, s = s/(GM), where s is a regularization scale 
with dimension of length; let us note that r\ 2 j s = r/s}. 
Although the logarithms do appear in terms that we 
are not able to regularize non-ambiguously by means of 
3-dimensional procedures, these procedures give unique 
values of the coefficients of the logarithms. The loga- 
rithms derived by us are in agreement with the results 
obtained earlier in [TJ [l2[ in the sense that by means of 
them we have recovered a computed earlier coefficient of 
logarithm in the formula for energy as a function of an- 
gular frequency along circular orbit [see Eq. (|3.3p below] . 
Instantaneous 4PN logarithmic terms were also derived 
in Ref. [26] from manifest non-local-in-time terms. 

We have computed but not displayed here the 4PN 
Hamiltonian in the general reference frame (which de- 
pends on some more undetermined coefficients) . We have 
checked that this Hamiltonian respects Poincare invari- 
ance in a way described in detail in Ref. [27]. Namely, 



we have constructed the most general template for the 
4PN center-of-mass vector (with all but one coefficients 
unspecified and with terms proportional to ln(ri2/s) 
present) and we have checked that it is possible to fit the 
coefficients of this vector in such a way that all Poincare 
algebra relations are fulfilled with 4PN accuracy. As a 
"by-product" of this checking we have uniquely deter- 
mined the quartic in momenta and proportional to v 3 
coefficient. 

The numerical coefficients not yet computed in the 
4PN center-of-mass Hamiltonian (|3.1j) . 12 different ones 
in total, need more complicated regularization proce- 
dures, likely even beyond pure dimensional regularization 
as successfully done e.g. in [2]. 

Making use of the Hamiltonian (|3.ip we have com- 
puted energy of binary system as a function of angular 
frequency along circular orbits (details of analogous com- 
putations performed at the 3PN order are given in Ref. 

28]). Along circular orbit n • p = and (p) 2 = j 2 /r 2 , 
where j = r x p (with j = |j|) is the reduced total an- 
gular momentum of the binary system. From equation 
dH(r, j 2 /r 2 )/dr = one gets the relation r — r(j) which, 
after substiution back to the Hamiltonian, defines energy 
E along circular orbits as a function of j, E = E(J). 
Then the angular frequency u> of circular orbits can be 
computed from the relation ui = (GM)~ 1 dE /dj . This 
relation can be perturbatively inverted to give j as a func- 
tion of uj and finally E = E(u). Instead of u one often 
uses the dimensionless PN parameter x = (GMw/c 3 ) 2 / 3 . 
We have also assumed that the regularization scale s can 
be related with the period P along circular orbits, s — cP 

12]. Making then use of the third Kepler's law at the 
Newtonian level one shows that s — 2nc~ 2 x~ 3 / 2 . Ap- 
plying the above described procedure to the Hamiltonian 
+ H ipn + ^2PN + ^3PN + -H4PN [the reduced Hamilto- 
nians from to H^pn can be found in Eq. (3.6) of [Uj], 
where one has to put w s t a tic = and Wkinctic = 41/24, see 
0], the 4PN-accurate binding energy of the system can 
be put in the form 



E(x; v) = - ^1 + eip N (^) x + e 2PN H x 2 + e 3PN (u) x 3 
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+ (e 4 PNM + ^In^) x 4 + O((v/c) 10 )^j , 



where the fractional corrections to the Newtonian energy at different PN orders read 
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(3.4a) 
(3.4b) 
(3.4c) 
(3.4d) 



In the above formula the two 4PN coefficients, at and 
v 4 , were computed for the first time. The 4PN coefficient 
c\ linear in v was computed numerically in Ref. 1291 . Its 
approximate value equals [see also Eq. (3.1) in j22| ] 



c x = 153.8803. 



(3.5) 
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